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Abstract 


The  "classical"  turbulence  approach  to  the  flow  of  air 
in  a  vegetative  canopy  is  reviewed.  Three  separate  "classical" 
theories  are  described.  A  comparison  is  made  between  experi¬ 
mentally  measured  canopy  flows  and  a  velocity  function  derived 
from  the  theory  of  structured  continua.  The  agreement  between 
theory  and  observations  is  good,  but  not  unexpected,  due  to  seven 
degrees  of  freedom  in  the  fitted  velocity  profile.  Proposals  for 
future  study  for  a  more  rigorous  test  of  the  theory  are  discussed. 
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I.  Introduction 

The  purpose  of  this  paper  is  to  apply  the  theory  of  structured 
fluids  to  the  turbulent  flow  of  air  in  a  vegetative  canopy.  The  model  ap¬ 
plied  here  is  the  flow  between  two  plates  of  fluid  containing  deformable 
structures.  It  is  believed  that  this  represents  the  first  attempt  to  com¬ 
pare  the  results  of  experiment  with  the  theory  of  structured  fluids. 

First,  the  classical  approach  to  canopy  flow  is  discussed. 
Then  the  theory  of  structured  fluids  is  reviewed  and  the  general  solu¬ 
tion  to  the  parallel  plate  problem  presented.  This  is  followed  by  a 
comparison  between  experimental  measurements  and  the  theory  of 
structured  fluids.  Finally,  the  results  are  summarized  and  discussed. 

II.  Classical  Approach  to  Canopy  Flow 

The  mean  flow  of  air  in  vegetative  canopy  has  been  studied 

’classically"  by  Ordway,  Ritter,  Spence  and  Tan  [l  J.  They  consider 

the  problem  essentially  the  same  as  the  classical  sliding  plate  problem: 

"At  the  upper  edge  of  the  flow,  say  z  z  we  see  that  the  outside  tur- 

h 

bulent  flow  effectively  drags  the  canopy  flow,  analogous  to  a  moving 

wall  with  velocity  U.  ..."(’ll. 

n 

The  basic  differential  equation  of  canopy  flow  is  derived  in 
[  1  ] .  The  principal  result  is: 
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whe  re 

z  =  vertical  coordinate, 

U  =  molecular  viscosity  coefficient  (r  . 

e(z)  =  eddy  viscos  ity  coefficient  (a  function  of  z h 

A (z )  =  total  vertical  plane  area  of  the  plants  comprising 
the  canopy  per  unit  volume  (a  function  of  z), 

I 

~  ''raP  coefficient  of  the  plants, 

p  mass  density'  oi  air, 

U  -  mean  horizontal  wind  velocity  (a  function  of  zl. 

The  physical  interpretation  of  this  equation  is  that  the  diver¬ 
gence  of  the  Reynold's  stress  is  equal  to  the  fluid  drag, 

l 

Before  equation  (1  I  can  be  evaluated,  e(z).  A ( r. > ,  and  must 

be  determined.  For  realistic  values  of  these  functions  equation  (1  1  can- 
no..  be  solved  analytically'  and  it  was  necessary  to  resort  to  a  numerical 
technique . 

Cionco  (2,  3 !  has  treated  this  problem  in  a  slightly  different 
manner.  Whereas  [l  )  was  eonrerned  with  the  behavior  of  a  ’transit*  r 
coefficient "  throughout  the  canopy  [2,  3j  concerned  itseli  w  ith  the  mix¬ 
ing  length  and  its  properties  throughout  the  canopy.  The  canopy  flow 
equation  rt:  rived  in  [  t!  is: 
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where 

f  =  mixing  length  of  the  canopy  (a  function  of  z), 

c 

X  =  z/h. 

C  and  G  are  two  parameters  introduced  by  Cionco,  their 
definitions  are  quite  lengthy  and  will  not  be  presented  here. 

Equation  (2)  may  be  interpreted  physically  in  the  same  man¬ 
ner  as  equation  (1  ).  It  must  also  be  evaluated  numerically.  The  pro¬ 
cedure  is  described  in  (3). 

A  purely  empirical  treatment  of  the  problem  was  given  by 
Uchijima  and  Wright  [4  j.  In  this  approach  the  authors  postulated  a 
velocity  ‘'unction: 

U  -  U,  +  V  In  ~  -  B ( z ) ,  0  z  h  (3) 

h  o  h 


w  he  r  e 


A  -  a  proportionality  constant  (a  Junction  of  U,  ), 
o  '  1  “ 

B  ( i- )  the  effect  of  the  plant  community  on  the  wind 
(a  function  of  •■), 
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The  authors  took  measurements  of  U  and  adjusted  B(z) 
so  that  the  equation  would  be  satisfied.  In  effect  B(z)  represents 
the  difference  between  the  measured  values  of  U  and  the  first  two 
terms  on  the  right  side  of  equation  (3),  It  would  appear  that  this 
procedure  allows  equation  (3)  to  agree  with  any  experimental  data. 

III.  Structured  Fluid  Approach 

Theory  of  structured  continua 

Basically  this  approach  treats  the  air  and  plants  as  a  struc¬ 
tured  continuum.  The  premise  introduced  here  is  that  the  bulk  fluid 
velocity  of  the  structured  fluids  approach  may  be  identified  with  the 
mean  velocity  of  the  classical  turbulence  approach. 

A  continuum  may  be  described  as  a  medium  in  which  the 
field  quantities  (displacements,  stresses,  velocities,  etc.  )  are  piece  - 
wise -continuous  functions  of  the  coordinates  of  the  material  points 
and  time. 

A  structured  continuum  is  a  continuum  that  contains  cer¬ 
tain  mathematical  abstractions  called  "structures".  At  present  there 
seems  to  be  a  number  of  possible  physical  interpretations  for  these 
structures,  see  [5,6,7]  for  example.  However,  none  of  these  have 
been  formalized.  Perhaps  some  formal  physical  interpretation  will 
evolve  in  the  near  future.  It  is  clear  as  pointed  out  by  Truesdell  and 
Toupin  [8]  that  with  or  without  such  an  interpretation  the  mathematics 
of  the  theory  retain  their  validity. 


The  axioms  upon  which  the  theory  is  based  are: 

(a)  Conservation  of  Mass, 

(b)  Balance  of  Momentum, 

fc)  Balance  of  Moment  of  Momentum, 

(d)  Conservation  of  Energy, 

(e)  Principle  of  Entropy. 

The  discussion  of  structured  fluids  presented  here  has 
been  necessarily  brief.  For  the  details  and  foundations  the  interested 
reader  is  referred  to  [5,6,758,9,  10,  1. 1 ,  12]. 

Constitutive  Equations  and  General  Theory 

When  forces  are  applied  to  any  material,  axioms  (a)  through 
(e  )  must  hold.  However,  when  the  same  system  of  forces  is  applied 
to  different  materials  the  materials  resp  nd  differently.  In  other 
words,  the  same  system  of  forces  applied  to  different  materials  qive 
rioe  to  different  strains  or  rates  of  strain.  The  relation  between 
stress  and  strain  is  referred  to  as  the  constitutive  equations. 

The  constitutive  equations  for  a  structured  fluid  have  been 
developed  in  [1  1],  and  the  parallel  plate  problem  has  been  solved  for 
a  fluid  containing  rigid  structures  [12]  and  for  a  fluid  containing  non- 
rigid  structures  [13],  The  formulation  leads  to  3  coupled  ordinary 


differencial  equations: 


f 


h 


(a+b)  D ' -t  2(B  -E?)  D  =  0, 

(4) 

(a-b)  V  ’  -  E  U1  4  2Z3  V  =  0, 

(5) 

tA2-E3)  U"  4-  2E2D'  +  2E  V'  -  P, 

(6) 

where 

U  -  mean  velocity  of  the  fluid, 

V  =  rate  of  rotr.tion  or  vorticity  of  the  structures, 

D  =  deformation  of  the  structures, 

P  =  the  pressure  gradient, 

a,b,  B^,  E^,  E^,  A.,  are  viscosity  coefficients. 

Equations  (4)  and  (5)  are  statements  of  axiom  (c )  and  equation  (6)  is 
a  statement  of  axiom  (b).  Equation  (4)  if  a  result  of  considering  de¬ 
formable  structures. 

4  ■ 

It  should  be  pointed  out  that  a  and  b  have  dimension  L  T 
■vhile  B^,  E^,  and  A^  have  dimension  L  T  .  Kirwan  f  1 3 j  dis¬ 
cusses  the  constraints  on  the  viscosity  coefficients  due  to  the  second 


i 


law  of  thermodynamics. 
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The  solutions  to  the  above  set  of  equations  with  z  nor¬ 
malized  were  derived  for  the  following  boundary  conditions  in  [l3] 

(a)  U(z  =  0)  =  U(z  -  1 )  r  0 

D  (z  =  0)  =  -D(z  =  1 )  =  Dq 
V(z  =  0)  =  -  V(z  =  1 )  =  VA 

(b)  P  =  0 

U (z  =  0)  =  0 

U (z  -  1)  =  u  >  0 
D  (z  =  0)  =  D(z  =  1  )  =  Dq 
V  (z  =  0)  =  V(z  =  1)  =  Vq 


The  solution  to  this  system  of  equations  with  general  boundary 

* 

conditions  for  the  lower  plate  fixed  and  the  upper  one  moving  is 


U  = 


..  ..  WV  .  , 

1  )  +  X  +  - r  (Lz 

z  0i„ 


2E  M 

+  — -  r(Kz) 

5 


(DQ  ;  Dt)  L (Lz) 


2E  N 

R(Kz), 

5 


(7) 


❖ 

The  vertical  coordinate,  z,  has  been  normalized  for  convenience. 
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V  =  Mf(Kz)  +  2Y jD  +  P^f  ~-l)  +  Nh(Kz)  +  ~  , 


(8) 


Here: 


D  ~  (Dq  "  Dj  )  f(Lz)  +  (DQ  +  D} )  h(Lz). 


u  A  =  Uft  =  0 
z-  0  0 


(9) 


U  ,  =  u. 

z- 1  1 


V  „  =  V 

z=  0  0 


V  ,  =  V, 

z=  1  1 


D*=o  =  Do  ■ 


D  =  D  , 
z~  1  1  ' 


L  =  2(B2  -  E2)/(a  +  b)  , 


K  =  -2E3A2/(a  +  b)  a5  , 

4  Y2 

x  =  ui  +  cT  T  S(L)  (Do  +  Di)  +  e2ns(k)  , 

5 


N  =  eIa5  ,V0+  Vl|-Yia5<D0  +  Dl,-“2J  -  2V1S(L)(D0  +  D1), 
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M  =  (V 


0 


vi>- VDo-Di>  + 


tf/vz1  _  [£.9sA.Y(1-z?  +  cosh  (yz)J 
1  ~  y  sinh  y 

S(va)=  [coah  (yz)  -  cosh  y(l-z)] 

' !  \  sinh  v  ’ 


r(\z)  =  g(-yz)  -  g(v)  , 


R(Yz)  =  S(yz)  +  S(y)  , 


f(yz)  =  [sinh  y(I-z^  -  sinh  (Yz)]/sinh  y  , 


h(yz)  =  [sinh  (yz)  +  sinh  y(1  -z)]/sinh  y  , 


Yj  =  E2K2/A2(L^  -  K2)  , 

Y2=2(E2  +  E3Y1)  * 

C  =  a5  +  2E3S(K)  , 


ll  is  orthwhile  noting  that  instead  of  reducing  the  problem 
to  3  coupled  ordinary  dif  "ential  equations,  one  equation  composed  of 
higher  order  derivatives  of  U  could  have  been  obtained. 
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IV.  C omparison  Between  Theory  and  Experiment 

The  first  step  in  the  comparison  between  experimental  data 

and  the  theory  of  structured  continua  was  to  determine  if  the  velocity 

function  equation  (7)  is  capable  of  producing  velocity  profiles  similar 
,  * 

to  tnose  measured  .  In  this  paper  only  the  velocity  function  U(z) 
will  be  treated,  A  short  cut  in  this  comparison  is  found  via  the  method 
of  Least  Squares.  The  triai  function  we  would  like  to  use  is: 

U  =  Pj(z^  -  1)  +  B^z  +  p^r  (Lz)  +  p^r (kz)  -i-  p^RfLz)  +  p^R(kz)  (10) 

The  primary  problem  is  to  determine  p  ,  p  p  ,  p  ,  p 

1  4  j  *r  j 

p,  ,  K  and  L  from  observations.  Although  for  each  profile  there  were 

o 

only  six  independent  measurements,  thirty -six  additional  points  aloi.g 
the  profile  were  interpolated  for  use  in  the  Least  Squares  procedure. 

The  determination  of  K  and  L  make  this  a  nonlinear  least  square  prob¬ 
lem.  While  computational  methods  are  in  existence  for  this  [l4,  15] 
a  quicker  and  easier  procedure  was  used.  Equation  (10)  was  treated 
as  a  linear  trial  function  with  values  for  K  and  L  specified.  Likely 
values  of  K  and  L  were  used.  Following  this  procedure  a  new  trial 
function  was  introduced  so  as  to  facilitate  programming: 


The  data  for  this  comparison  was  supplied  by  Lemon  and  Allen  of  the 
New  York  State  College  of  Agriculture,  Cornell  University. 


1 1 


+  P4e 


-Lz 


.  a  -Kz  •  Kz 
+  P,  e  +  P.,  e 


(1 


It  should  be  noted  that  a  seventh  degree  of  freedom  has  been  added. 

This  is  because  we  are  no  longer  constraining  U(z  =  0)  to  be  0  due 
to  a  lack  of  data  in  the  lower  portions  of  the  canopy.  Tests  indicate 
that  this  constraint  has  little  or  no  effect  on  the  results.  The  results 
of  these  computations  are  illustrated  graphically  in  Figures  1  through 
9.  The  solid  lines  are  the  measured  profiles  and  the  broken  lines  are 
the  fitted  profiles.  In  the  cases  where  only  a  solid  line  is  plotted  the 
fitted  curve  was  so  close  to  the  measured  curve  that  no  distinction 
could  be  made  graphically.  A  total  of  approximately  twenty  thousand 
runs  were  made.  The  correlation  coefficients  in  every  case  were  be¬ 
tween  .  990  and  .999. 

As  can  be  seen  the  agreement  between  equation  (7)  and  the 
data  is  quite  good.  However,  it  was  found  that  just  as  good  as  fit  could 
be  obtained  with  a  sixth  degree  polynomial.  It  was  also  found  that  the 
fit  with  L  =  0  (rigid  structures)  and  a  fourtn  order  polynomial  were 
nearly  as  good. 


The  fact  that  an  arbitrary  trial  function  with  the  same  number 
of  degrees  of  freedom  as  the  velocity  function  derived  from  the  theory 
of  structured  continua  can  fit  the  data  so  well  is  to  be  expected  since 
the  degrees  of  freedom  of  the  trial  function  is  slightly  greater  than  the 
number  of  independent  data  points.  The  advantage  of  using  the  velocity 
function  derived  here  is  that  it  is  derived  from  a  physical  theory. 

V.  Discussion  of  Results 

As  can  be  seen  from  Figures  1  through  9,  the  velocity  func¬ 
tion  derived  via  the  theory  of  structured  continua  can  quite  accurately 
r  produce  experimentally  measured  mean  velocity  profiles  in  various 
vegetative  canopies.  However,  as  mentioned  previously,  the  fact  that 
the  trial  function  used  in  the  least  squares  routine  has  seven  degrees  of 
freedom  pretty  much  assures  a  good  fit  to  the  data.  It  is  worth  noting 
that  neither  the  log  wind  nor  the  parabolic  trial  function  could  reproduce 
the  experimentally  measured  profiles. 

lo  summarize  there  are  three  methods  of  reproducing  experi¬ 
mentally  measured  'Canopy  Flows": 

(a)  "Classical  Approach", 

fb)  Theory  of  Structured  Continua, 


(c1  Arbitrary  trial  function  with  enough  degrees  of  freedom. 
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Of  the  three,  the  last  seems  least  attractive  because  it  is 
completely  arbitrary  and  can  be  related  to  no  physical  concepts.  The 
first  one  suffers  from  what  might  be  called  "creeping  empiricism". 

Such  features  as  a  height  dependent  eddy  viscosity  and  the  complica¬ 
tion  of  an  unwieldy  differential  equation  that  cannot  be  solved  analyti¬ 
cally  make  this  method  unattractive.  The  remaining  method  is  the 
theory  of  structured  continua,  There  are  several  advantages  to  this 
approach. 

(a)  It  yields  a  velocity  function  that  is  relatively  easy  to  evaluate. 

(b)  There  are  rone  of  the  height  dependent  parameters  to  be  deter¬ 
mined. 

(c)  This  seems  to  be  a  logical  extension  of  the  continum  hvpo- 
pothesis  in  the  sense  that  the  structures  are  considered  to 
be  continuums  within  continuums. 

To  be  sure,  the  theory  of  structured  continua  involves  difficulties. 

For  one  thing,  there  is  no  physical  description  or  feeling  for  what  the 
structures  are.  Without  this  description  the  deformation  and  voriticity 
fields  cannot  be  measured  and  hence  only  one -third  of  the  theory  a  n  be 


investigated  experimentally. 
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VI.  Recommendations  for  Future  Research 

It  is  apparent  that  the  observations  used  in  this  study  are  not 
detailed  enough  to  adequately  test  the  structured  fluid  approach.  In  order 
to  obtain  more  conclusive  results  two  major  problems  must  be  solved  first 

(1)  better,  more  detailed  flow  measurements  must  be  otained; 

(2)  some  formal  physical  explanation  for  the  "si  uctures"  must 
be  developed. 

With  the  solutions  to  the  above  two  problems  in  hand  one  could 
use  the  differential  equations  (4,5,  6)  to  determine  all  of  the  viscosity 
coefficients  that  permeate  the  theory. 

At  this  point,  a  clearly  written  orderly  summary  of  the  theory 
of  structured  fluids  would  be  quite  beneficial.  The  abstract  formalizations 
have  reached  the  point  where  they  must  be  stated  concisely  and  put  in 
their  proper  perspective  so  that  the  interested  researcher  will  have  a 
jumping  off  point  from  which  he  can  wade  through  the  material  already 
in  print  on  the  subject. 
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Type  of  Canopy 

zH(cm) 

Lth (cm/sec)  JK 

L 

Remarks 

Japanese  Larch 

1040 

149 

3.0 

3,5 

Can't  distinguish 
graphically  between 
data  and  fit. 

Japanese  Larch 

1040 

348 

3.5 

3.0 

Can't  distinguish 
graphically  between 
data  and  fit. 

Oats 

125 

61 

4.5 

3.0 

Can't  distinguish 
graphically  between 
data  and  fit. 

Oats 

125 

169 

3.0 

3.5 

Slight  distinction 

Soy  Beans 

150 

98 

3.0 

3.5 

No  distinction 

Soy  Beans 

150 

246 

2.5 

3.0 

Slight  distinction 

Corn 

300 

81 

4.5 

3.0 

Slight  distinction 

Corn 

300 

300 

3.5 

3.0 

No  distinction 

Sunflowers 

225 

225 

3.5 

3.0 

No  distinction 

Sunflowers 


225 


225 
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